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Построен гамильтониан кристалла, испытывающего структурный переход мартенситного типа, в котором система конфигура- 
ционных возбуждений описывается как двухуровневая квантовая система (квантовая система псевдоспинов). 


1. Введение 

Экспериментально обнаружено, что кристалл, 
неустойчивый относительно структурного перехо- 
да исходная - конечная структура, вызванного как 
изменением температуры, так и внешней силы, в 
окрестности структурного перехода находится в 
предпереходном состоянии [1, 2]. При описании 
сегнетоэлектрических, сегнетоэластических, бес- 
порядок-порядок переходов широко используются 
модельные представления о двухямном кристалли- 
ческом потенциале, псевдоспине, успешно разра- 
батываются теоретические методы исследования 
таких систем [3-6]. Однако для описания сегнето- 
электрических переходов второго рода [3] при раз- 
ложении потенциальной энергии атомов достаточ- 
но ограничится двухчастичным взаимодействием 
псевдоспинов, в то время как для случая мартенси- 
тных превращений первого рода необходимо учи- 
тывать трехчастичное и четырехчастичное взаимо- 
действие псевдоспинов. 

Для теоретического описания структурного 
превращения мартенситного типа и предпереход- 
ного состояния сформулирована микроскопиче- 
ская модель [7, 8], описывающая конфигурацион- 
ные возбуждения структурнонеустойчивого кри- 
сталла как квантовую систему псевдоспинов. Усло- 
вие возникновения предпереходного состояния в 
кристалле заключается в том, что внешнее воздей- 
ствие (изменение температуры, механическая си- 
ла), стимулирующее структурный переход исход- 
ная - конечная структура, уменьшает площадь гор- 
ба, разделяющего минимумы двухямного кристал- 
лического потенциала атома. Это приводит к суще- 
ственному увеличению квантового туннелирова- 
ния атома (при низких температурах), тепловых 
перескоков (при высоких температурах) и умень- 
шению асимметрии двухямного кристаллического 
потенциала. Возникает неустойчивость состояния 
исходной кристаллической решетки с асимметрич- 
ным двухямным потенциалом относительно воз- 
никновения предпереходного состояния решетки с 
симметричным двухямным потенциалом. 

В этой работе построен гамильтониан кристал- 
ла, испытывающего структурный переход мартен- 
ситного типа, в котором система конфигурацион- 
ных возбуждений описывается как двухуровневая 
квантовая система (квантовая система псевдоспи- 
нов) [7, 8]. 


2. Псевдоспиновый формализм и гамильтониан 

структурнонеустойчивого кристалла 

Из полученных в [7, 8] оценок следует, что если 
узел исходной структуры находится от сопряженно- 
го узла конечной структуры на расстоянии порядка 
6<10~ 9 см (меньше амплитуды нулевых колебаний 
атома), а площадь горба, разделяющего левый и пра- 
вый минимумы потенциала, менее К 2 6=1,5-10~ 22 
эрг-см, то необходимо учитывать квантовое тунне- 
лирование атома между сопряженными узлами ис- 
ходной и конечной структур. То есть наряду с малы- 
ми колебательными смещениями атома внутри ле- 
вой потенциальной ямы (фононами) в двухямном 
потенциале появляются дополнительные квантовые 
смещения атомов (туннелирование) в определен- 
ном направлении и на определенное расстояние - 
дискретные конфигурационные степени свободы. 
Следовательно, волновая функция атома должна за- 
висеть не только от непрерывной пространственной 
координаты х, но и от одной дискретной перемен- 
ной, указывающей значение проекции псевдоспина 
на некоторое выбранное направление в простран- 
стве псевдоспина, например, ось?. Для нашего слу- 
чая двухуровневой системы волновая функция ато- 
ма будет иметь вид спинора ЧДх,^), который пред- 
ставляет собой совокупность двух различных функ- 
ций координат четной Т(х,+1/2)=Т + (х) и нечетной 
Т(х,— 1/2) =ѵ РДх), отвечающих различным значе- 
ниям ^-компоненты псевдоспина. Оператор псев- 
доспина при применении его к волновой функции 
ЧЧх,5 г ) действует только на переменную 5\ Для опе- 
раторов псевдоспина выполняются обычные ком- 
мутационные соотношения 


где а,[і,у=х,у,і в пространстве псевдоспина, а,Ъ - 
нумеруют ячейки, - подрешетки кристалла в ко- 
ординатном пространстве. Для каждой двухуровне- 
вой системы любой оператор частиц (эрмитова ма- 
трица второго порядка) может быть выражен через 
операторы псевдоспина 1/2 
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В представлении четной и нечетной волновых 
функций Т+, ЧД, одночастичная потенциальная 
энергия атома в ячейке а, подрешетке у в симме- 
тричном двухямном потенциале будет иметь вид [3]: 

Ѵ\КІ) = Псо 81 р 

где Нео - расщепление энергий четного и нечетно- 
го состояний атома в двухямном потенциале, К] - 
радиус вектор координаты атома в ячейке. 

Далее следует учесть один важный аспект рас- 
сматриваемой задачи, в основе которого лежит экс- 
периментальный факт - структурное превращение 
в кристалле. Этот факт можно интерпретировать 
так, что в исходной структуре атом кристалла нахо- 
дится в левой яме двухямного потенциала, а в ко- 
нечной структуре в правой яме двухямного потен- 
циала. Следовательно, в структурнонеустойчивом 
кристалле в качестве основного состояния двуху- 
ровневой системы псевдоспинов можно выбрать не 
дублет из четной и нечетной волновых функций Т + , 
Т_, а дублет из волновых функций <р ъ ер к , локализо- 
ванных в левом и правом положении потенциала, 
соответственно. Очевидно, что собственные функ- 
ции ер г , ер к являются симметричной и антисимме- 
тричной комбинацией волновых функций ЧД, Ч* 

Ѵь =(1/2) 1/2 (Ч / + + Ч / _), 

<Р Я =( І/2) 1/2 (Ч' + -Ч'_). 

В представлении локализованных волновых 
функций ер ь ср к одночастичная потенциальная энер- 
гия атома в ячейке а, подрешетке у в симметричном 
двухямном потенциале будет иметь вид [3]: 

= ( 1 ) 

Физический смысл псевдоспиновых операто- 
ров 8-, 8- проясняется в представлении локализо- 
ванных состояний (р,, ср к . Оператор 8- ( I ) характе- 
ризует разницу заселенностей симметричного и 
антисимметричного состояний, а оператор 8- ха- 
рактеризует разницу между заселенностями левого 
и правого положений двухямного потенциала. 

При построении гамильтониана структурноне- 
устойчивой решетки необходимо учесть, что свя- 
занные с переходом особенности превращений 
мартенситного типа проявляются в их термодина- 
мике и динамических свойствах не на очень боль- 
ших частотах, меньших или порядка характерных 
фононных (например, дебаевской частоты кри- 
сталла а>дяа10 13 с~‘) и много меньших характерных 
частот электронной подсистемы. Поэтому при 
описании свойств структурного превращения мар- 
тенситного типа можно исходить из обычного для 
кристаллов адиабатического приближения [3, 9], 
считая, что состояние кристаллической решетки 
полностью определяется заданием координат и 
скоростей атомов, и не рассматривать явно элек- 
тронных степеней свободы. Тогда гамильтониан 
кристаллической решетки будет иметь вид [3,9]: 


Я = Х(Р^/2т ; )+К(КІ), (2) 

а) 

Ѵ аі - оператор импульса в ячейке а, подрешетке у, 
Й(К') - потенциальная энергия атомов. Суммиро- 
вание идет по всем атомам решетки (координатам 
ячеек а и номерам подрешеток у). 

В окрестности структурного превращения мар- 
тенситного типа атом, наряду с малыми колеба- 
ниями относительно положения равновесия (фо- 
нонами) в одной из ям двухямного потенциала, мо- 
жет время от времени как совершать тепловые пе- 
рескоки, так и туннелировать через потенциаль- 
ный барьер в другую потенциальную яму. Частота 
такого перескока может быть меньше или оказать- 
ся сравнимой с дебаевской частотой кристалла. 
Поскольку равновесные свойства системы не зави- 
сят от времен релаксации, необходимо учитывать 
переходы атома между ямами потенциала. 

Изменения двухямного потенциала, вызванные 
малыми колебаниями атома относительно положе- 
ния равновесия (фононами), будут небольшими. 
Поэтому в статическом приближении можно пре- 
небречь малыми колебаниями и при вычислении 
термодинамических свойств решетки рассматри- 
вать только равновесные конфигурации атомов, а 
влияние фононов учесть по теории возмущений. 

Для симметричного двухямного потенциала с 
двумя положениями равновесия два возможных 
значения равновесной координаты атома К^іѵ+Ь, 
аналогично [3] можно записать с помощью опера- 
тора Паули 8' аі в виде: 

( 3 ) 

где величина х аі - радиус вектор средней координаты 
атома в ячейке а , подрешетке у исходной структуры. 
Далее, следуя [3], предположим, что потенциальную 
энергию атомов К(К') (2) можно записать в виде сум- 
мы одно-, двух-, трех- и четырехчастичных взаимо- 
действий: 

П(К^) = Х^(Щ)+(1/2)Х^'(С^) + 

а,і а,Ъ,і,і 

+(1/3) X К*(К',К',К*) + 

а,Ь,с,і,і,к 

+(1/4) X ^ и (кі,щ,к*,к;). (4) 

а,Ъ,с4,і,},к ,т 

Используя метод [3], для равновесной конфигу- 
рации атомов, определяемой значениями 8% (3), 
потенциальную энергию (4) можно записать, огра- 
ничиваясь членами до четвертого порядка по 8 % , в 
виде: 

ГЮ = І>^ -(1/2) I - 

а,] а,Ъ,і,І 

-(і/з) х / # (г л-чте- 

а,Ъ,с,і,І,к 

-(1/4) X (5) 
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где константа X выражается через потенциал 
Х' аа: = Р(г ш +аЬ,,г 4у +а Ъ у ) с а,а-± 1, константа У* - 
через потенциал ^^ а .= РХ г оі+«Ь, 5 г і; +аЪу,г Л +а'\) с 
а, а', а" =+1, а константа / № ” - через потенциал 
•ЛЬ ѵ= ^(г о; + аЬ;,т Ь] + а Ъ р г ск + а "\>„г іт + а ”Ь И ) с 
а,а',а",а"-±\ следующим образом: 

^(г а „г,) = (1/4 )<УІ + ѴЦ - VI - VI), (6) 

^(г а ,,г 4я г Л ) = (1/8)х 

х(К^ +к^ -Ѵ®_ -Ѵ + І к + +ѵ;і -ѵ т + +ѵ * ), (7) 

Кікт Кх ЬР ѵл т )= 


= (1/16)(К + ^ + -Г + ^_ - V І кт + +ѴІ ікт _ +ѵ*% - 

ту у кт I уукт . ту у кт ту у кт . ту і/ кт , ту у кт 

-++- -+-+ ‘ -+ — "++++ '"+++- ''++-+ 

уукт уукт уукт уукт уукт ^ 


( 8 ) 


Из (5) видно, что включение двух-, трех- и че- 
тырехчастичного взаимодействия между псевдос- 
пинами вызывает асимметрию двухямного потен- 
циала, то есть возникает кристаллическая решетка 
с исходной структурой. 

Учтем малые колебания атомов около положе- 
ния равновесия и' (фононы) для завершения по- 
строения полного гамильтониана (2). Для этого по- 
тенциальную энергию К(К)) (5) разложим по степе- 
ням смещений и', имеющим вид: 


При структурных превращениях мартенситного 
типа меняется как размер, так и форма кристалла. 
Эти изменения принято описывать введением бес- 
конечно малых величин - компонент тензора 
упругих деформаций Поэтому выражения 
(10-12) можно переписать в более удобном виде: 

Н рИ =\сІѵ{^[Р 2 а (г)/2р] + (\/2)Х аРуд е^(ѵ)а'' в {г)}, (13) 

а 

н м = I КУ/, + + а /2)5^х; \ау<* . ( 14) 

Здесь р - средняя плотность кристалла, Р а - 
компонента плотности импульса поля малых упру- 
гих деформаций (фононов), Х ік1т - тензор модулей 
упругости кристалла [9], а р - расстояние между 
атомами по Д - оси координат, е“Р - тензор упругих 
малых деформаций кристалла в месте расположе- 
ния атома в ячейке а, подрешетке у, связанный с 
компонентами его малых смещений и" соотноше- 
нием [9] 


^ =(1/2 )[ди:/дх р +д и ^/дх а ]« 
~ ( 1 / 2 ) [ \У“ д ГІ + И// ] ехр ( і(\ г аі ) , 


< = = (9) 

Из формул (9) и (4) видно, что замена К-'^К^+и/ 
соответствующая учету малых колебаний в (4), эк- 
вивалентна формальной замене г а] — >г а/ +и ; а в (5). Раз- 
лагая правую часть (5) с точностью до членов перво- 
го и второго порядка по и' [3], получим гамильтони- 
ан структурнонеустойчивой решетки: 

н = Н рк + V (К° у )+ Н^, (10) 

я *=Е< р */2*,)+ 

+(1/2) X (1 / 2)< ; («: - < №-и»\ (11) 

а,Ы,і 

Н ім = I К^ + ^)5"-(1/2)5да“1(«; - и*). (12) 

а,Ъ,і,і 

Здесь Н рИ - гамильтониан малых колебаний ре- 
шетки (фононов), К(К/) - гамильтониан взаимо- 
действия псевдоспинов, Н м - гамильтониан псев- 
доспин-фононного взаимодействия, которое изме- 
няет как квантовое туннелирование атома, так и 
асимметрию двухямного потенциала. По повто- 
ряющимся индексам а, /3, у, нумерующим декарто- 
вы компоненты векторов, в выражениях (10-12) и 
ниже подразумевается суммирование от 1 до 3. Ве- 
личины А“[‘, В?, У" выражаются через производные 
потенциалов (6-8) следующим образом: 

= д/дх а {В\Г}, {А*? } = д 2 /д*?дх р {А*}, 
АЧг аі ,Гу) = (1/4 )(УІ + ѴЦ + VI + Ѵі), 
ВЧг ш ,г Ь] ) = (1/4 )(УІ - ѴЦ + VI - ѴЦ 


Интегрирование в (13) ведется по всему объему 
кристалла. Здесь, аналогично [3, 6] применяется 
полуклассический подход, при котором динамика 
псевдоспинов описывается квантовомеханически, 
а поле малых упругих деформаций (фононов) - 
классическим образом. По этой причине гамильто- 
ниан (13), в отличие от ( 14), представляет собой не 
оператор, а числовой функционал. 

Исследование структурных переходов, описы- 
ваемых гамильтонианом (10), является достаточно 
сложным, поэтому при рассмотрении статических 
и динамических свойств превращений мартенси- 
тного типа будем использовать метод самосогласо- 
ванного поля [3]. В этом методе поле, действующее 
на атом со стороны остальных, заменяется на сред- 
нее, определяемое из условия самосогласования. 
Поскольку структурный переход заключается в по- 
явлении отличного от нуля среднего смещения, 
равного й / (0,5-<5|>), из положения равновесия в 
левой яме в положение равновесия правой ямы 
двухямного потенциала, то для построения после- 
довательных приближений метода самосогласо- 
ванного поля удобно выделить в операторе псев- 
доспина 5 ’е=< 5'|>+(5|— <5|>) среднее значение 
псевдоспина <3%> и флуктуацию В нуле- 

вом приближении метода самосогласованного по- 
ля, называемом приближением молекулярного по- 
ля, взаимодействием флуктуаций пренебрегают 
[3], и гамильтониан (10) можно заменить на эф- 
фективный: 

н м = Нрк + Ѵ м (кУ)+Н м іп1 , (15) 
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=-2Х'§» + 

+(1/2)Х Л(г т .,г 4 .)<5^><5^> + 

+(2/3) X І *(\, ѴЪХЧ х5* >< Х/ > + 

+(3/4) X ЯПь.Ь.Ъ.Тм) < « х -5? х « >< •%, >. 

(16) 

ь - = (-й©-йп;,о,-йп* + 

+ Х / ‘'( г -* г *)<‘ х *> + 

+ X / “( г -> і ѵ> г л) < ^ >< > + 

+ X < 5* х 5Д >< 51 », (17) 

Ъ,с,й,],к,т 

ш \і = X ( ^ + С / 2) < ^ ^ )а р ^ , 

Ъ,і 

= (18) 

Ъ,і 

Среднее значение псевдоспина в элементарной 
ячейке й, подрешетке / определяется как [3]: 

(8 Ш > = 8р[8 ш ехр-(/1// м )]/8р[ехр-()37/ м )]. 

Очевидно, что описание структурнонеустойчивой 
решетки при помощи гамильтониана (15-18) по- 
прежнему остается сложным (для полиморфного, а не 
изоморфного мартенситных превращений), посколь- 
ку в (15—18) присутствуют члены, описывающие 
взаимодействие конечных конфигурационных сме- 
щений из различных подрешеток. Однако поскольку 
структурный переход - экспериментальный факт, 
причиной которого является такое изменение взаи- 
модействия псевдоспинов, которое приводит к впол- 
не определенному изменению структуры, то в рамках 
предложенной модели это можно описать, наложив 
ограничения на константы взаимодействия: 

•П^Л,) » Л(г аі ,г ы), /“(г ш .,г и ) » /“(г ш .,г 4 Д 

В п(?аП Г Ъ,)» Щ( Г аП*Ъі) ДЛЯ І * У, 

/"Чг,, , Т ы , Г а . ) » Р к (г ш . , г 4 . ,т л ) для ІФІФк, 

К ‘ Ш ( Г аП Г Ы’ Г С Л)» 

для іФ]ФкФт. (19) 

Условие (19) означает, что одинаковые конфи- 
гурационные смещения атомов сильнее взаимо- 
действуют между собой, чем различные конфигу- 
рационные смещения атомов, в результате чего эф- 
фективный гамильтониан (16-18) распадается на 
независимые гамильтонианы подрешеток: 

Ѵ м ( Ц ° а І )+ Н м ы = ~Х I* А + 

+(1/2)Х/"(г,,г 4 ,.)<^.><^.> + 

+(2/3) X 1“Р а М,) < 5* >< 5 ; х > + 

+(3/4) X К т Р„г ы ,т сЛ ) < % х % >< 52 х 52 >, (20) 

а,Ь,с4,і 


К =(-йю-Ш:,.,0,-Ш:,. + Х/“(г а ,.,г 4 ,.)< 52> + 

ь 

т«)<^х^> + 

Ь,с 

+Х ^‘“Чг ш ,г 4 ,,г„г л )< 55 X X х 52 >), (21) 

Ь,с,(1 

ШХ. = Х(Д“+(1/2)<55>/“)а^, 

Ь 

( 22 ) 

ь 

Выполняя термодинамическое усреднение 
(20-22), получим среднюю энергию системы псев- 
доспинов: 

<Н М + ]<$*> + 

+[Х (5“ + (1/ 2) < 55 > Г й )а,4* ] < 52 > +Н р „ + 
-(1/2)ХЛ г „,г к )<5‘х^>- 

а,Ьі 

-(1/3) X ^(1^,0 < 52 х 5* х 52 > - 

а,Ь,с,і 

-(1/4) X 7Г'(г в ,., г и , г с ,., г л ) < 52 х 52 х 52 х 52 > . 

а,Ъ,с,<Ц 

Благодаря аддитивности гамильтониана (20-22) 
при изучении особенностей статистики и динамики 
структурнонеустойчивого кристалла можно ис- 
пользовать упрощенную модель, в которой в эл- 
ементарной ячейке находится только один атом ре- 
шетки, а суммирование по / можно опустить. 

4. Обсуждение результатов 

Предложенная модель позволяет заключить, 
что систему конфигурационных возбуждений в 
структурнонеустойчивом кристалле необходимо 
описывать как квантовую систему псевдоспинов 
(квантовую двухуровневую систему). В окрестно- 
сти структурного перехода исходная - конечная 
структура внешнее воздействие уменьшает пло- 
щадь горба, разделяющего минимумы двухямного 
потенциала атома. Это приводит к возникновению 
эффектов квантового туннелирования, теплового 
перескока атома и уменьшению асимметрии двух- 
ямного потенциала, что открывает возможность 
переходов из узлов исходной решетки в узлы ко- 
нечной решетки (конфигурационных смещений) и 
возникновения предпереходного состояния. Под 
предпереходным состоянием кристалла понимает- 
ся такое конденсированное состояние кристалла, в 
котором атом решетки, вследствие эффектов кван- 
тового туннелирования и теплового перескока, 
полностью делокализован в симметричном двух- 
ямном потенциале, то есть когда вероятность об- 
наружить атом в узле исходной и конечной струк- 
туры одинакова. Таким образом, по мере увеличе- 
ния внешнего воздействия ангармонические эф- 
фекты нарастают, а конфигурационные смещения 
атомов из узлов решетки увеличиваются и начина- 
ют взаимодействовать между собой. В результате 
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этого происходит потеря устойчивости решетки в 
определенных кристаллографических направле- 
ниях, возникает ближний порядок смещений (ста- 
тические смещения) атомов, протекает структур- 
ное превращение. Оценим порядок величин 

/=2>.,г 4 ,г е) , 

Ь Ь,с 

К = 1 1 К ( Т а’ Г Ь’ Г с’ Г Л 
Ь,с,сі 

а 11 В" =а^В"(г а ,г ь ), а» Г =<Р (г а ,г ь ). 
ь ь 

Поскольку величина конечных смещений ато- 
мов Ь при превращениях мартенситного типа мала 
по сравнению с межатомным расстоянием 
|х"-х 4 “| =а а , то потенциалы 

Ь Ъ,с 

Ь,с4 

можно разложить по степеням Ь. Тогда из (5-8) вид- 
но, что К имеет четвертый порядок, I - третий, /, 
а р ] а - второй, а р Б" - первый, 4, Д, А 2 - нулевой по- 
рядок степени Ь, т.е. степени отношения 
(й/а^ЗТО^.-.З-КТ 1 ). Если считать, что А 0 , А,, Л 2 име- 
ют обычный атомный порядок величины (ед. эВ), то 

а 1 ' В" =(3-10- 2 ...3-10“ 1 )4, 

/,вѴ“ =(9 10 4 ...9 10 2 )Д, 

/ = (27-10“ 6 ...27-10“ 3 )Д, 

^ = (8Ы0Л..8Ы0- 4 )Д. 

Подставляя значения этих величин в среднюю 
энергию системы псевдоспинов и считая, что к 
кристаллу приложена внешняя сила, стимулирую- 
щая структурное превращение {(^ву^уо, 

|^(б" +(1/2) < 8* > /“ )а р е аР ]> 0) и вызывающая 

Ь 
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статическую деформацию ^«З-КМ.-З-Ю -1 , полу- 
чим изменения знака как у двух-, так и у трехча- 
стичного взаимодействия псевдоспинов. Это сви- 
детельствует о том, что внешняя сила способна вы- 
звать превращения мартенситного типа путем пере- 
хода исходная структура - предпереходное состоя- 
ние - конечная структура, наблюдаемые экспери- 
ментально [1,2]. Именно по причине, что внешнее 
воздействие (изменение температуры, внешняя си- 
ла) может изменять знак двух- и трехчастичного 
взаимодействия псевдоспинов в окрестности пре- 
вращения мартенситного типа, при разложении по- 
тенциальной энергии нельзя ограничится двухча- 
стичным взаимодействием, как в случае сегнетоэ- 
лектрического [3], и необходимо учитывать трех- и 
четырехчастичное взаимодействие псевдоспинов. 

Разумно предположить, что эффект квантового 
туннелирования (тепловых перескоков) атомов яв- 
ляется физической причиной когерентного пове- 
дения кристалла при низко(высоко)температурном 
структурном переходе мартенситного типа, стиму- 
лированным внешним воздействием. Вызванная 
этими эффектами полная делокализация волновой 
функции атома позволяет рассматривать предпере- 
ходное состояние в кристалле, как суперпозицию 
исходной и конечной структур. 

Анализ полученных результатов, теоретических 
расчетов [7, 8] и экспериментальных данных [1,2] 
позволяет сделать вывод, что предпереходное со- 
стояние кристаллической решетки, испытывающей 
мартенситный структурный переход, при изменении 
как температуры [1], так и внешней силы [2], харак- 
теризуется спонтанным возникновением конечных 
(конфигурационных) смещений атомов и сопровож- 
дается упругой деформацией, обладает низкой 
устойчивостью. Поэтому слабоустойчивое состояние 
деформируемой кристаллической решетки [2], мож- 
но считать аналогом предпереходного состояния 
кристалла, испытывающего мартенситный переход. 
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